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We analyze the equal-time Bethe-Salpeter quark wave functions of the pion in various models. We 
discuss how the quenched lattice QCD results with delocalized pion interpolators can be identified 
with the coarse grained wave functions, typical of low-energy effective models. Actually, we find 
that one-loop chiral quark models predict that pseudoscalar and tensor wave functions have the 
same shape, while the axial component is more extended. These facts are accurately confirmed by 
the lattice. We also show how the transversity information, relevant for the light-cone physics, can 
be straightforwardly obtained from the equal-time rest-frame lattice calculations. This remarkable 
relation provides a way to extract, for instance, the equal-time holographic wave functions and 
compare them, quite favorably, to the lattice calculations. 
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1. Introduction 

Hadronic wave functions encode important information on bound states in strong-interaction 
physics; in particular, they provide the amplitude for a composite hadron to have quarks in a 
given momentum state or, equivalently, at a certain space-time distance. Heavy quarks obey non- 
relativistic quantum mechanics and conserve particle number. While our understanding and intu- 
ition is based on wave functions, as a matter of principle the wave functions cannot be directly 
measured experimentally. One must instead resort to form factors, decay widths, or momentum 
distributions. Moreover, for light quark systems particle creation may occur, demanding a field- 
theoretic framework where further complications arise. Relativistic invariance requires that one 
uses the conventional Bethe-Salpeter (BS) amplitudes with a fixed number of the quark field oper- 
ators, a reminiscent of the approximated parton picture point of view, emphasized by the light-cone 
approaches [|J]. Color gauge invariance requires additional inclusion of the link operators 

For the pion, the spontaneously broken chiral symmetry is a basic dynamical ingredient in 
the determination of its nonperturbative quark structure. It appears via the pertinent axial Ward- 
Takahashi identities pi]. These important constraints are implemented in relativistic field-theoretic 
chiral quark models, such as the Nambu-Jona-Lasinio (NJL) model (for a review see, e.g., [Q]). 
The regularization, introducing the physical cut-off, needs to be carefully handled not to spoil the 
relativistic, gauge, and chiral symmetries. 

On the other hand, lattice QCD solves the bound state problem in a fundamental way. It is 
thus possible to make a first-principle nonperturbative determination of the wave functions, but 
at the expense of breaking the continuum symmetries, such as the Lorentz invariance and, quite 
often, chiral symmetry, due to the finite lattice spacing. The axial Ward-Takahashi identities can be 
exactly implemented on the discrete Euclidean lattice as shown by Ginsparg and Wilson [§] (see 
Ref. Q for a recent practical implementation), enabling realistically small pion masses. 

In the present contribution we show the analysis of the pion wave functions from the quenched 
lattice QCD ^ and make the comparison to various hadronic models. In spite of the very dissim- 
ilar appearance and nature of these approaches, we will provide the conditions under which this 
comparison may be undertaken. We also address in more detail the light-cone issues with the help 
of the transversity relations. 

2. Bethe-Salpeter Amplitudes 

The BS vertex or wave function of the pion is given by 



where q{x) are spinor field operators carrying flavor and color, and \n a (p)) is the pion state with 
the Cartesian isospin index a and the on-shell four-momentum p, p 2 = ra|. While chiral quark 
model calculations are naturally formulated in the momentum space, the basic objects in Euclidean 
lattice calculations are the point-to-point correlation functions. These quantities are gauge and 
renormalization-group invariant at all Euclidean times, which basically correspond to off-shell pro- 
cesses. At large Euclidean times only the on-shell states contribute to the correlation functions, as 




(2.1) 
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well as the off-shell violations of the gauge invariance disappear. Inverting the Fourier transforma- 
tion we get (in the isospin limit of m u = mj) 

(0\T {q(x)q(0)} \n a (q)} = ^-k^p + Wa + ic^q^T], (2.2) 

J% 

where the wave functions a = P,A,T, depend on the Lorentz-invariant variables x 2 , x-q, and 
q 2 = m\. The quantities *F a are vertex functions in the BS equation, and as such are finite and 
undergo x-independent multiplicative renormalization. Thus, the ratios *¥ a (x) (0) become cut- 
off independent, as the cut-off is removed, which on the lattice means a — > 0. 



The definition of Eq. (2.2) is completely satisfactory for chiral quark models. In QCD, how- 
ever, it is only gauge-invariant in the fixed-point Fock-Schwinger gauge, x^A^{x) = 0, where the 
standard derivatives, d 11 , and the covariant derivatives, = d 11 +igA^, coincide. On the lattice 
the gauge fixing has the problem of the Gribov copies, as there is no complete gauge fixing. On the 
other hand, Elizur's theorem prevents non-vanishing vacuum expectation values of gauge variant 
operators in the physical Fock space. 

Non-gauge invariant operators can be made gauge invariant by joining them with a link oper- 
ator, however, as a result the path-dependence sets in. Furthermore, gluons carry momentum in the 
pion and different gauge-invariant definitions yield different results (see Ref. [§] for a discussion 
on various possibilities). For definiteness, we choose a straight-line path and undertake a smearing 
procedure. This derealization improves the signal-to-noise ratio for the measured hadron corre- 
lators, as the interpolating operators have a larger overlap with the desired state. Local operators, 
in contrast, do not take into account the spatial extension of the hadron. The usefulness of the 
smearing process lies also in the fact that the overall thickness of the flux tube in the probe is con- 
trolled by the number of the smearing steps. In addition, the method is computationally simple. 
The resulting fat link also reduces the high-energy fluctuations and the path dependence, such that 
we deal with a coarse-grained wave function. The procedure naturally finds its counterpart in the 
low-energy effective chiral quark models. In a previous work the quenched lattice calculations 
of the pion have been worked out along these lines. The quenched approximation contains all the 
leading -A^, and hence the qq, Fock state components. Thus we expect that quenched calculations 
describe the large-A c motivated models 1 . 

3. Transversity relations 

The relativistically invariant BS amplitude has the representation 

(0\T{q(x)q(0)}Mq)) = iWa C dae^ 2a -^ x 

Jo 

x [-^p(a,x 2 ) + $> A (a,x 2 ) -Ha^q^xy^Tia^ 2 )], (3.1) 



where a is the Feynman parameter. As a matter of principle, all scalars such as *P fl in Eq. ( J2.2[ ) de- 
pend on the kinematic variables x 2 , x ■ q, and q 2 , thus we are free to choose any form of kinematics. 



'The quenched lattice calculations also contain a piece subleading in N c , which is actually suppressed for heavy 
quarks; pion loops are l/Af c -suppressed, although not all of the l/iV c -contributions originate from the pion loops WM. 
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In the rest-frame kinematics we have xo = and (qo, q) = (m n ,0), whence x = — x , x ■ q = 0, and 
q 2 = m\. On the other hand, in the infinite-momentum-frame kinematics (go, q) = (\J m \ + p 2 ,p z ), 
with p z — > oo. Thus on the light-cone surface, xq = z, one has x 2 = — x|, and x ■ q — > 0. 

The Lorentz invariance allows us to relate the rest-frame calculation to the transverse-coordinate 



dependence in the light-cone wave functions. Simply, by comparing Eqs. ( |2.2D and ( p.ip we find 

,-\ „ 

x ¥ a {x 2 ,x.q)= / doW a {a,x?-)e- iq - x V- a - l) . (3.2) 
Jo 

We may identify the Feynman parameter a with the Bjorken x- variable, ^Bj = OC. Then, for the 
chosen kinematics q + x~ = q x = 0, we have *P ET (— r 2 ,0) = J dxBj^^ixBj, — r 2 ), where ET and 
LC denote the equal-time, and light-cone wave functions, respectively. In the argument of *F ET 
one takes r 2 = x 2 , the distance squared, while in the argument of *P LC we need to use r 2 = x|, the 
transverse distance squared. Therefore 

¥f (-r 2 ,0) = jf'dxBj ^(xbj -4)L r=r - (3-3) 

That way the connection between the ET and LC wave functions has been established 2 . 

Although generating the autonomous connection between LC and ET wave functions seems at 



first glance hopeless [10], a similar trans versity relation has recently been deduced for scalar parti- 



cles. Another transversity property for the Generalized Parton Distributions (GPDs) was suggested 



for the nucleon [ |1 1| | and the pion [|12|], allowing a frame-independent definition of probability (un- 
like the more conventional Breit-frame definition). In the case of the pion it reads 

P LC (b) = J'dxq(x,h) = J ^e^ b F v (-ql)=P ET (r)\ r=hMbl (3.4) 

where q(x, b) is the off-forward diagonal GPD and Fy (/) is the pion form factor. In the phenomeno- 
logically successful Vector Meson Dominance (VMD), where Fy(t) = M v / (My — t), one gets 

P ET (r) = M v K (rM v )/(2n) ~ e*/^ . (3.5) 

Besides these relations, it would also be useful to verify the ET-LC transversity connection directly 



on the lattice. While there exist transverse lattice calculations ni3n (see [14] for a review), their 



focus is placed on the Distribution Amplitude, *P(a,0) = <p(ot), leaving out the x± dependence. 
4. Chiral Quark Models vs Lattice 



We evaluate the correlation function of Eq. ( \1.2\ > in a chiral quark model (for a review see, 
e.g., |0]). Disregarding for the moment the regularization, an instructive way to determine the 
pion wave function in a chiral quark model is by exploiting the axial Ward-Takahashi identity. 
It relates the quark propagator, S(p), and the vertex function corresponding to the axial current, 
J^ ,a (x) = ^q(x)y^ Y5T a q(x) , with the irreducible vertex r^' a (/? + <?,/?): 



S(p + q)- l Y5\?a + Y5\?aS{pT l = q ^> a (p + q , p) . (4.1) 



2 Note that we keep only the qq components of the BS amplitude, thus we do not account for the possible emergence 
of the higher Fock-state components, pertinent to the dynamical nature of the boost. 
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Figure 1: The components of the rest-frame equal-time pion wave function, normalized to unity at the 
origin, x V a {r)/ x i' a (Q), for S, A, and T channels, evaluated in the NJL model at M = g nqq f„ = 300 MeV, 
and compared to the quenched lattice data rf7h, The displayed points of the lattice data are at m K = 345 and 
575 MeV, while the NJL model calculation includes also the case of the physical pion mass. 



In the NJL model, the spontaneous breaking of the chiral symmetry generates a constituent quark 
mass, M, given by the so-called gap equation. As the result, S(p) = l/(p — M), such that 



„ q^2M 

t-—~r 



(4.2) 



The pole at q 2 = indicates the Goldstone boson nature of the pion. The pion wave function is 
extracted from the pion pole as an unamputated vertex function, 



l + q-M \f n 



M 



l-M 



(4.3) 



where the Goldberger-Treiman relation at the quark level, g nqq =M/f K , can be read off. With the 
Feynman trick, the result becomes particularly simple in the chiral limit q 2 = m\ — > 0, yielding 

y P (a,x 2 ) = M[-2(2a- l)d# +5% +M 2 ] d M2 A(M,x), 
^ A (a,x 2 ) = M 2 d M iA(M,x), 
■$V(a,x 2 ) = -2d x 2d M 2A(M,x), 



where we have introduced the free scalar propagator in the coordinate space, 
A(M,x) -- 



d 4 p e'P* MK^MV^x 2 ) 



(2tt) 4 p 2 -M 2 

From the previous formulas we find (for m n = 0) the relations 



\n 2 \J-x 2 



e- Mr 
7372"' 

e -Mr 



(4.4) 



(4.5) 



(4.6) 



where ^0 and ^1 are the modified Bessel functions and "reg" means a regulator. The asymptotic 
behavior at r — ^ 00 is independent of the regulator and implies a longer tail in the A channel than 
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Figure 2: Axial equal-time holographic wave functions, normalized to unity at the origin, for the quark 
masses m = 330 MeV (dashed black) and m = 4 MeV (solid red), obtained via the transversity relations and 
compared to the lattice calculation with m K = 575 MeV (boxes) and 345 MeV (circles) [pp. 



in the P and T channels. The exponential fall-off of the pion wave functions is consistent, up to 
a power in r, with the probabilistic estimate in Eq. ( |3~3| ) with My = 2M, M being the constituent 
mass. At short distances as A(x 2 ,M 2 ) ~ 1 jx 2 the results are divergent demanding regularization. 

The NJL model with the (twice subtracted) Pauli-Villars (PV) regularization applied to an 
observable A amounts to the replacement M 2 — > M 2 + A 2 , followed by the subtraction 



A| rcg - A ( A 2 -(),- .,l(A-)+A 2 ^S. (4.7) 



Results of the calculation, shown already in [f7|], are given in Fig. ffl. As we see, the agreement is 
excellent. At the farthest lattice point, r = 0.6 fm, the ratio rW^ (r) /*J/p(r) approaches a constant. 

The quantitative agreement of the NJL model with the data is not trivial. For instance, in 
the Spectral Quark Model (SQM) [ 15] only qualitative matching is achieved. In this approach the 
regularization is introduced by replacing the mass M by the "spectral" mass CO, using g nqq = co/f K 
and integrating over CO with a suitable weight p(ffl), which depends on the vector and scalar meson 
masses (My and M$), and also on a specified contour in the complex co-plane. That way, for 
instance, VMD of the pion form factor can be built in with M v = 24 7i 2 f 2 /N c . For m n = we get 
the results 

«F„(r) _ ¥ r (r) _ = e -M v r,2 + (4 . 8) 



¥ P (0) ¥r(0) ' *a(0) 

Again, ^ is more extended than *J/p and Wj, due to the presence of an extra power in r. In the case 
of SQM we would get good fits of the lattice data for My jl = 505(30), 520(20), 530(14) MeV for 
the subsequent values of m n = 345, 475, 575 MeV. A simple quadratic extrapolation in m n to the 
chiral limit yields M v /2 = 493(20) MeV, a too high value as compared to My = m p = 770 MeV. 



5. Holographic wave functions vs lattice 



The transversity relations (3.3) can be used to deduce the ET wave functions from the LC wave 
functions. As we have shown, ET wave functions can be computed on Euclidean lattices upon a 
suitable coarse graining of the gauge link operator. As an example, we take here the holographic 
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wave functions inspired by the AdS/CFT correspondence and recently brought in connection with 
the LC wave functions [16] as a first approximation to QCD (for a review see e.g. [17]). The 
basic idea was to relate the QCD-LC Hamiltonian by using the scaling relation between the two- 
dimensional vectors £ = y / x(l — x)b (valid for massless quarks) and assume that the interaction 
also depends on this scaling variable. Further elaborated models with finite quark masses are 
introduced, assuming the replacement k 2 — > + m 2 [18]. However, by doing so it is not obvious 
whether m corresponds to the current or constituent mass of the quark. For that reason, masses m = 
4 MeV and m = 330 MeV are explored. A full discussion of models is carried out in Jl9|], where 
the soft-wall with a positive dilaton background seems phenomenologically preferred. Actually, 
besides the good quality of the mass spectrum, M 2 LS = 2no(n+L+S/2) with 4fc 2 = 2no , the pion 
arises as a massless mode (corresponding to n = L = S = 0). Unlike the chiral models, this is not 
linked to the spontaneous breaking of the chiral symmetry, apparently not manifest in the light-cone 
dynamics. Thus, we identify the BS axial component with the holographic wave function [16, 18] 



—=y/x(l -x) exp 
\ % 



K 2 , 



-x)b\ 



ii r 



2k 2 x(1 -x) 



(5.1) 



Calculations can be undertaken analytically for massless quarks. Following [ |18| ] we fix the pion 
weak decay, 7t + — > }X + V, and the neutral pion decay, 71° — > 2y, from the conditions 



1 



dx^>\ c {b L ,x) 



2^ Jo v-^lbx=o.- 2V 3 

respectively, with f n = 92 A MeV. This yields A 2 = 2/n 2 = 0.2 and K = An \ 
Using the transversity relation the ET (holographic) wave function reads (in the massless case, 
0) 



Ak 
T6 



2^/k [ l dxd 2 b ± *&f(b ± ,x) 
Jo 



4An z 



yj> = 

2/3 f n =950 MeV. 



m 



ET/ 



^T(O) 



,-2/^ 2 r 2 /3 



[/ (2/^V/3)-/ 1 (2/^ 2 r 2 /3)] =1 



(5.2) 



where I n (z) are the modified Bessel functions. Note that upon the use of the relation My = 
24n 2 f%/N c the small-r expansion reproduces the SQM result, Eq. (4.8), exactly. The asymptotic 
behavior at large distances has the form 



(5.3) 



^T(O) 8/ 3 7T 7 /2 ^ ^ 

which displays the exponential fall-off, similarly to the chiral quark models, Eq. ( |4.6[ ) for NJL and 
Eq. ( pT^ ) for SQM, however the powers of r are different, exhibiting different dynamics in the 
models. The parameters A and k for m = 4 MeV and m = 330 MeV in Eq. ( |5Tl ) are taken as in 
[18]). As we can see from Fig. ||, the agreement is quite good and the lattice data hardly allow to 
discriminate between the mass values except at r = 0.6 fm. Of course, it would be very useful to 
pin-down the correct long distance behavior from the lattice above 0.6 fm. The calculation of the 
other (higher-twist) components within the holographic approach would also be highly desirable. 

Our numerical study shows that the lattice data in Fig. || can be best fitted with m ~ 100 MeV, 
similar to the current quark mass used in the NJL calculation at high values of m % . In particular, 
for m n = 575 and 345 MeV we have used m = 140 and 51 MeV, respectively, in rough agreement 
with the Gell-Mann-Oakes-Renner relation —m(qq) = f„m„. Thus, the interpretation of m in the 
holographic models as the current quark mass seems consistent with this comparison. 
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6. Conclusions 

The presented calculations show that the Euclidean lattices can be successfully used to coarse 
grain the wave function over short-distance scales, where the gluon degrees of freedom are inte- 
grated out. A direct comparison to wave functions from various hadronic models not only becomes 
meaningful, but in some cases very successful. A quite unexpected result concerns the utility of 
our calculations to determine the transversity information (the dependence on the transverse co- 
ordinates), relevant for the light-cone physics; the infinite-momentum frame pion wave functions, 
integrated over the Bjorken-x, coincide in the impact-parameter space with the equal-time rest- 
frame wave functions. This relation provides a way of checking the wave functions for models 
genuinely formulated in the LC variables. As an example, we have carried out this analysis for 
holographic models. 

We thank Sasa Prelovsek and Luka Santelj for their collaboration in the lattice calculation [pi]. 



References 



[1] S.J. Brodsky, H.-C. Pauli and S. S. Pinsky Phys. Rept. 301 (1998) 299-486 Qhep-ph/9705477| ]. 

[2] H. Suura Phys. Rev. D17 (1978) 469. 

[3] H. Pagels and S. Stokar Phys. Rev. D20 (1979) 2947. 



[4] E. Ruiz Arriola Acta Phys. Polon. B33 (2002) 4443-4479 [ |hep-ph/0210007[ |. 
[5] P. H. Ginspai-g and K. G. Wilson Phys. Rev. D25 (1982) 2649. 



[6] C. GattringerP/zys. Rev. D63 (2001) 114501 Qhep-lat/0003005p . 

[7] W. Broniowski, S. Prelovsek, L. Santelj and E. R. Arriola Phys. Lett. B686 (2010) 313-318 



[0911.47051 



[8] J. W. Negele hep-lat/000702 6 



[9] T. D. Cohen and D. B. Leinweber Comments Nucl. Part. Phys. 21 (1993) 137-150 



[hep-ph/9212225] 



[10] G. A. Miller and B. C. Tiburzi Phys. Rev. C81 (2010) 035201 [0911 .3691], 



[11] G. A. Miller Phys. Rev. Lett. 99 (2007) 112001 [0705.2409] 



[12] W. Broniowski, E. R. Arriola and K. Golec-Biernat Phys. Rev. D77 (2008) 034023 [0712 .1012], 



[13] S. Dalley and B. van de Sande Phys. Rev. D67 (2003) 1 14507 [hep-ph/0212086] 



[14] M. Burkardt and S. Dalley Prog. Part. Nucl. Phys. 48 (2002) 317-362 [hep-ph/0112007]. 



[15] E. Ruiz Arriola and W. Broniowski Phys. Rev. D67 (2003) 074021 [hep-ph/0301202). 



[16] G. F. de Teramond and S. J. Brodsky Phys. Rev. Lett. 102 (2009) 081601 JO 8 9 . 4 8 9 9p . 
[17] S.J. Brodsky and G. F. de Teramond 



1009 . 4232 



[18] A. Vega, I. Schmidt, T. Branz, T. Gutsche and V. E. Lyubovitskij Phys. Rev. D80 (2009) 055014 



[0906. 1220| 



[19] S. S. Afonin 1001 .3105 



8 



